Identification of a rare SARS-CoV-2 XL hybrid variant in wastewater and the subsequent discovery of two infected individuals in Nevada.

The identification of novel SARS-CoV-2 variants can predict new patterns of COVID-19 community transmission and lead to the deployment of public health resources. However, increased access to at-home antigen tests and reduced free PCR tests have recently led to data gaps for the surveillance of evolving SARS-CoV-2 variants. To overcome such limitations, we asked whether wastewater surveillance could be leveraged to detect rare variants circulating in a community before local detection in human cases. Here, we performed whole genome sequencing (WGS) of SARS-CoV-2 from a wastewater treatment plant serving Las Vegas, Nevada in April 2022. Using metrics that exceeded 100× depth at a coverage of >90 % of the viral genome, we identified a variant profile similar to the XL recombinant lineage containing 26 mutations found in BA.1 and BA.2 and three private mutations. Prompted by the discovery of this rare lineage in wastewater, we analyzed clinical COVID-19 sequencing data from Southern Nevada and identified two cases infected with the XL lineage. Taken together, our data highlight how wastewater genome sequencing data can be used to discover rare SARS-CoV-2 lineages in a community and complement local public health surveillance.

Use of wastewater surveillance for early detection of Alpha and Epsilon SARS-CoV-2 variants of concern and estimation of overall COVID-19 infection burden.

A decline in diagnostic testing for SARS-CoV-2 is expected to delay the tracking of COVID-19 variants of concern and interest in the United States. We hypothesize that wastewater surveillance programs provide an effective alternative for detecting emerging variants and assessing COVID-19 incidence, particularly when clinical surveillance is limited. Here, we analyzed SARS-CoV-2 RNA in wastewater from eight locations across Southern Nevada between March 2020 and April 2021. Trends in SARS-CoV-2 RNA concentrations (ranging from 4.3 log10 gc/L to 8.7 log10 gc/L) matched trends in confirmed COVID-19 incidence, but wastewater surveillance also highlighted several limitations with the clinical data. Amplicon-based whole genome sequencing (WGS) of 86 wastewater samples identified the B.1.1.7 (Alpha) and B.1.429 (Epsilon) lineages in December 2020, but clinical sequencing failed to identify the variants until January 2021, thereby demonstrating that 'pooled' wastewater samples can sometimes expedite variant detection. Also, by calibrating fecal shedding (11.4 log10 gc/infection) and wastewater surveillance data to reported seroprevalence, we estimate that ~38% of individuals in Southern Nevada had been infected by SARS-CoV-2 as of April 2021, which is significantly higher than the 10% of individuals confirmed through clinical testing. Sewershed-specific ascertainment ratios (i.e., X-fold infection undercounts) ranged from 1.0 to 7.7, potentially due to demographic differences. Our data underscore the growing application of wastewater surveillance in not only the identification and quantification of infectious agents, but also the detection of variants of concern that may be missed when diagnostic testing is limited or unavailable.

Genomic surveillance of Nevada patients revealed prevalence of unique SARS-CoV-2 variants bearing mutations in the RdRp gene.

Patients with signs of COVID-19 were tested through diagnostic RT-PCR for SARS-CoV-2 using RNA extracted from the nasopharyngeal/nasal swabs. To determine the variants of SARS-CoV-2 circulating in the state of Nevada, specimens from 200 COVID-19 patients were sequenced through our robust sequencing platform, which enabled sequencing of SARS-CoV-2 from specimens with even very low viral loads, without the need of culture-based amplification. High genome coverage allowed the identification of single and multi-nucleotide variants in SARS-CoV-2 in the community and their phylogenetic relationships with other variants present during the same period of the outbreak. We report the occurrence of a novel mutation at 323aa (314aa of orf1b) of nsp12 (RNA-dependent RNA polymerase) changed to phenylalanine (F) from proline (P), in the first reported isolate of SARS-CoV-2, Wuhan-Hu-1. This 323F variant was present at a very high frequency in Northern Nevada. Structural modeling determined this mutation in the interface domain, which is important for the association of accessory proteins required for the polymerase. In conclusion, we report the introduction of specific SARS-CoV-2 variants at very high frequency in distinct geographic locations, which is important for understanding the evolution and circulation of SARS-CoV-2 variants of public health importance, while it circulates in humans.